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prolific pulsed source of electron and muon neutrinos as well as muon antineutrinos. 
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INTRODUCTION 



The question of detecting and exploiting neutrinos from both terrestrial and extra terrestrial 
sources has become central to physics and astrophysics. In recent years, in particular, exploiting the 
neutral current interaction for detecting neutrino induced nuclear recoils has become fashionable. 
Indeed due to the coherence of all neutrons in the nucleus one gets very large cross sections, which 
are proportional to the square of the neutron number (the coherent neutral current contribution due 
to protons is negligible). To exploit this feature, however, one should design large detectors with 
very low energy thresholds. Two possibilities have been considered: 

• Spherical TPC detectors filled with a noble gas for supernova neutrino detection. Such detec- 
tors can give us information about the neutrino source, since the neutral current interaction 
is not sensitive to detailed neutrino properties like neutrino masses, mixing etc, which appear 
in neutrino oscillations. It has also been pointed out that such detectors are robust, cheap 
and easy to maintain for years. So they can surely detect supernovae at distances of tens of 
kiloparsecs 

• The standard direct dark matter detectors aim at detecting nuclear recoils due to WIMP's 
(weakly interacting massive particles). Since the expected event rates are extremely low, one 
may have to worry about the rare background recoils due to neutrino sources, like the solar 
8 B neutrinos Q, Q. 

It is thus important to study the behavior of such detectors using well known terrestrial neutrino 
sources. One such source is Spallation Neutron Source (SNS) at the Oak Ridge National Laboratory, 
Oak Ridge, Tennessee, which is the most powerful source of pulsed intermediate-energy neutrinos 

The primary parameters of the SNS are as follows: 
Proton beam power on target: 1.4 MW 
Proton beam kinetic energy: 1.0 GeV. 
Average beam current on target: 1.4 mA 
Protons per pulse on target: 1.5 x 10 14 
Pulse repetition rate: 60 Hz 
7T+ per proton: 0.068 
7r + per second 6.14 x 10 14 
R(y e ) = Rfa) = 6.14 X 10 14 
The total neutrino flux at a distance r is: 

$(i/ e ,r)=$(i/^,r) = *(i/ M ,r) = — 5 (1) 



At a typical distance of 50 m we find: 



$(i/ e ,r) = $(i^,r) = $(^,r) = 1.95 x lOVn-V 1 
The normalized neutrino spectra can be described very well by the shapes exhibited in Fig. [1] 
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FIG. 1: The neutrino spectrum from stopped pions. The normalized dashed and solid curves correspond to 
v e and £ M respectively. Also shown is the discreet spectrum (thick vertical line). 



ELASTIC NEUTRINO NUCLEON SCATTERING 



The cross section for elastic neutrino nucleon scattering has extensively been studied. It has been 
shown that at low energies it can be simplified and be cast in the form: @,0| : 
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where tun is the nucleon mass, T/v its energy and gv, 9a are the weak coupling constants. Neglecting 
their dependence on the momentum transfer to the nucleon they take the form: 



gv = — 2 sin 
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In the above expressions for the axial current the renormalization in going from the quark to the 
nucleon level was taken into account. For antineutrinos g A — * —gA- To set the scale we write: 
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The nucleon energy depends on the neutrino energy and the scattering angle, the angle between the 
direction of the recoiling particle and that of the incident neutrino. In the laboratory frame it is 
given by: 



(7) 



(forward scattering). For sufficiently small neutrino energies, the last equation can be simplified as 
follows: 
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The above formula can be generalized to any target and can be written in dimensionless form as 
follows: 

2C0S 2 # Trecoil E v 

V = 7TTT"/ — \2 2^ ' y = > x » = ^ 

(1 + l/x v y - cos z 9 m recoi i m recoi i 

In the present calculation we will treat x v and y as dynamical variables, in line with CDM recoils. 
One, of course, equally well could have chosen x v and 9 as relevant variables. 
The maximum energy occurs when 9 — 0, i.e.: 

2 

VmaX ~ (1 + 1/^)2-1' (9) 

in agreement with Eq. (2.5) of ref. 0]. One can invert Eq. [5] and get the neutrino energy associated 
with a given recoil energy and scattering angle. One finds 
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The minimum neutrino energy for a given recoil energy is given by: 



, 0<6»<tt/2 (10) 
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in agreement with Eq. (4.2) of ref. 0. The last equation is useful in obtaining the differential cross 
section (with respect to the recoil energy) after folding with the neutrino spectrum 

COHERENT NEUTRINO NUCLEUS SCATTERING 

From the above expressions we see that the vector current contribution, which may lead to coher- 
ence, is negligible in the case of the protons. Thus the coherent contribution [8( may come from the 
neutrons and is expected to be proportional to the square of the neutron number. The neutrino- 
nucleus scattering can be derived in analogous fashion, ft can also be obtained from the amplitude 
of the neutrino nucleon scattering by employing the appropriate kinematics, i.e. those involving the 
elastically scattered nucleus and the substitution 



P P La , P 2 E A 



with q the nucleon momentum, A is the nuclear mass number and p the nuclear momentum. Under 
the above assumptions the neutrino-nucleus cross section takes the form: 



,2/1 T A,2 , /nr2 A/T 2\AmNT A , 



+ (M v - Ma) (1 - -fY + (Mi - My) £ " ] (12) 



where My and Ma are the nuclear matrix elements associated with the vector and the axial currents 
respectively and Ta is the energy of the recoiling nucleus. The axial current contribution vanishes 



for + =>■ + transitions. Anyway, it is negligible in the case of coherent scattering on neutrons. 
Thus Eq. (HH) is reduced to: 

(£r) =^P^ {NVm oh {T A ,E v ), (13) 

\ dT Aj wea k 2lT 

with 

F coh (T A ,E„) = F\f ) (l + (1 - T -±f - (14) 

where N is the neutron number and F(q 2 ) — F{T\ + 2ArriNTA) is the nuclear form factor. The 
effect of the nuclear form factor depends on the target, since the maximum recoil energy depends 
on the target (see Fig. [2|). 
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FIG. 2: The square of the nuclear form factor, F 2 (Ta), as a function of the recoil energy for A=131 (top) 
and A=40 (bottom). Note that the maximum recoil energy is different for each target. 



QUENCHING FACTORS AND ENERGY THRESHOLDS 

. The above results refer to an ideal detector operating down to zero energy threshold. For a real 
detector, however, as we have already mentioned, the nuclear recoil events are quenched, especially 



at low energies. The quenching factor for a given detector is the ratio of the signal height for a 
recoil track divided by that of an electron signal with the same energy. We should not forget that 
the signal heights depend on the velocity and how the signals are extracted experimentally. The 
actual quenching factors must be determined experimentally for each target. In the case of Nal the 
quenching factor is 0.05, while for Ge and Si it is 0.2-0.3. For our purposes it is adequate, to multiply 
the energy scale by an recoil energy dependent quenching factor, Q/ ac (T J 4) adequately described by 
the Lidhard theory 0]. More specifically in our estimate of Qu(Ta) we assumed a quenching factor 
of the following empirical form (Toj : 



Qf ac (T A ) = n 
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n ~ 0.256 , r 2 ~ 0.153 



(15) 



The quenching factors very much depend on the detector type. The quenching factor, exhibited in 
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FIG. 3: The quenching factor as a function of the recoil energy in the case of A=131 (top) and A=40 
(bottom). 



Fig. Ofor recoil energies of 131 Xe and Ar, were obtained assuming the same quenching of the form 
of Eq. (|15p . In the presence of the quenching factor as given by Eq.f fTB")) the measured recoil energy 
is typically reduced by factors of about 3, when compared with the electron energy. In other words 
a threshold energy of electrons of 1 keV becomes 3 keV for nuclear recoils. 



RESULTS 



For a heavy target, 131 Xe 

The neutrino-nucleus cross sections are shown in Figs|4]and[5l The quenching factor does not affect 
the rate (see Figs |4] and [5|) but it shifts the threshold downwards. In other words the unavailable 
phase space is on the left of the thin line with no quenching and the thick line with quenching. 
Clearly the effect of quenching is going to be larger when the recoil energy is smaller. Integrating 
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FIG. 4: The differetial neutrino nucleus cross section as a function of the recoil energy in the case of A=131 
for the continuous v e spectrum (top) and the continuous spectrum (bottom). The solid (dashed) curves 
correspond to no form factor (form factor) respectively. The fine vertical line corresponds to a threshold of 
3 keV. The phase space on its left is unavailable. Due to the presence of quenching the excluded phase space 
becomes larger (on the left of the thick line) . 

the differential rates shown Figs 2] and [5] at zerow threshold we obtain the total rates shown in table 
HJ The effect of a non zero threshold is exhibited in Figs [6] and [7J It is clear that both the effect 
of the inclusion of the form factor and of quenching are important, especially in the case of high 
energy threshold, but the effect of quenching is most important of the two. 
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FIG. 5: The same as in Fig. [4] for the discreet spectrum. 



TABLE I: The neutrino-nucleus total cross section in units of 10 41 cm 2 due to the neutral current 



target 


(no FF) 


(FF) 


(no FF) 


(FF) 


(no FF)) 


(FF) 


m Xe 


2.71 


1.98 


3.60 


2.46 


2.79 


1.79 


4U Ar 


0.221 


0.190 


0.294 


0.244 


0.178 


0.162 



For a light target, Ar 

The differential cross sections are now shown in Figs [5] and [5] Integrating the differential rates 
shown Figs [5] and [5] at zero threshold we obtain the total rates shown in table ITT1 The effect of a non 
zero threshold is exhibited in Figs [10] and [TTJ The effect of quenching now is less dramatic since 
the recoil energy for a light target is higher. The total number of events for a target of mass m in 
time t is given by 

R = 3.156 x 10 7 ^— L)a(v, (A, N))^-N A (16) 
lkg ly A t 

where Na — 6.02 x 10 23 is Avogadro's number, A t the mass number of the target, a(y, (A, N)) is 
the coherent neutrino-nucleus cross section and < &(^, L) is the neutrino flux (neutrinos per second 
per cm 2 ) in distance L from the source. Using 

§{v,L = 50m) = 1.96 x 10 6 cm- 2 s _1 

for each flavor we obtain the data of table HT1 



EXPERIMENTAL ISSUES WITH THE SPHERICAL TPC 



The SNS facility provides an excellent opportunity to employ and test the spherical TPC gaseous 
detector llj of volume V under pressure P and temperature T filled with a noble gas. This detector 




FIG. 6: The event rates at a given energy threshold divided by the corresponding ones at zero threshold 
as a function of the thrshold energy in the case of A=131 target for the continuous v £ spectrum (top) and 
the continuous spectrum (bottom). The solid (dashed) curves correspond to no form factor (form factor) 
respectively. The two lowest curves take into account to the quenching factor. 



has good resolution and very low energy threshold. In this case we obtain: 

t PV 
R = 3.156 x 10 7 — §{y,L)a(v, (A,N))—s(L,x) 
ly kT 

where the parameter s(L, x) depends on the shape of the vessel and the distance L of its center from 
the source. Using 

$(z/,i = 50m) = 1.96 x 10 6 cirr 2 s _1 

for each flavor we obtain the data of table Hill The function s(L,x) for a sphere of radius R with its 
center at a distance L from the source is given by 

T 2 i-L+n, r n i 



< R > R ^ - J^R^j L _ R ^ V + r 2 + 2 ,Lcosfl (i7) 




FIG. 7: The same as in Fig. [6] for the discreet spectrum. 



TABLE II: The event rates for nuclear recoils per ton per year due to neutral current neutrino-nucleus 
scattering for a detector at a distance L=50m away from the neutrino source with zero energy threshold. 



target 


(no FF) 


(FF) 


(no FF) 


(FF) 


(no FF)) 


(FF) 


all v 
(no FF)) 


all v 
(FF) 


m Xe 


7.70 


5.64 


10.3 


6.99 


6.21 


5.09 


24.0 


17.69 


4l, Ar 


2.06 


1.77 


2.74 


2.27 


1.66 


1.51 


6.46 


5.5 



The ratio s(R, R/L)/s(10, 0.2 is shown and Fig. [T2]. The experimental set up is shown in Fig. [T3] 
Using an appropriate field corrector the field inside the detector is spherically symmetric (see Fig. 

M). 

DISCUSSION AND CONCLUSIONS 

The total number of events per ton-y at 50 m is: 

i?(Xc) = 17.8(events per ton-y), R(Ar) — 5.5(events per ton-y) (18) 
One notices the following: 

1. The content of tables [LTllIIII can be understood by comparing it to other neutrino sources. 
The flux at 50 m is $ = 3<l>p e — 5.85 x 10 6 cm~ 2 s _1 . This is comparable to that of solar 
boron neutrinos — 5.0 x 10 6 cm~ 2 s _1 , but with a higher cross section, since the latter is 
characterized by lower average energy. Thus the relatively low flux and energy of the boron 
solar neutrinos cannot cause any background problems to dark matter searches 0. Anyway 
the number of neutrinos in a year is N v = 1.8 x 10 14 cm~ 2 . The corresponding quantities for 
supernovae neutrinos at a distance D — 3.1 x 10 22 cm, are [l[ 3?^ = 1.5 x 10 12 cm~ 2 , = 
1.0 x 10 12 cm~ 2 and $ Va = 6.3 x 10 n cm~ 2 for all other flavors a, i.e a total of 3.1 x 10 12 cm~ 2 . 
In other words, since the energy spectra are not very different, the event rate here is expected 
to be higher than that of a typical supernova in our galaxy [l|. 



2. The event rate per unit target mass predicted does not favor the heavy target according to 




20 40 60 80 100 120 140 

Ta -> keV 




20 40 60 80 100 120 140 

Ta -> keV 



FIG. 8: The differential neutrino nucleus cross section as a function of the recoil energy in the case of A=40 
for the continuous v e spectrum (top) and the continuous spectrum (bottom). The solid (dashed) curves 
correspond to no form factor (form factor) respectively. The fine vertical line corresponds to a threshold of 
3 keV. The phase space on its left is unavailable. Due to the presence of quenching the excluded phase space 
becomes larger (on the left of the thick line) . 



the expected N 2 rule implied by the differential cross section (compare Figs Q][S] and [SKI ) for 
the following reasons: 

• There fewer atoms in a kg of a heavy target. So one expects: 

N 2 

Rot—. 
A 

• The phase space available favors a light target, since, for a given neutrino energy, the 
nuclear recoil energy decreases with A. This kinematical advantage maybe partly offset 
by the fact that the effect of form factor becomes more important, if the recoil energy 
becomes higher. 

• Threshold effects may somewhat complicate the issue (compare Figs [6][7] and fl^lfTTT) . 
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FIG. 9: The same as in Fig. [8] for the discreet spectrum. 



TABLE III: The event rates for nuclear recoils per year due to neutral current neutrino-nucleus scattering 
for a detector with the shape of a sphere of radius R = 10 m with its center at a distance L = 50 m 
from the source(s(10,0.2)=1.008). The vessel is filled with gas under pressure P = 10 atm and temperature 
T = 300 °K. The detector is assumed to operate at zero threshold. 



target 


v R 
(no FF) 


(FF) 


(no FF) 


(FF) 


(no FF)) 


(FF) 


all v 
(no FF)) 


all v 
(FF) 


13i Xe 


1705 


1249 


2280 


1548 


1393 


1120 


5379 


3917 


4U Ar 


139.3 


119.7 


185.3 


153.5 


112.2 


102.1 


436.86 


375.3 



This ./V 2 dependence clearly applies whenever the amount of the target material for a given 
volume is determined by the pressure (see table Mil) in agreement with our earlier work on 
supernova neutrino detectors [l| . 
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FIG. 13: A sketch of spherical TPC detector to be employed in studying the coherent scattering of neutrionos 
eployed in this work. 
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FIG. 14: The Equipotential surfaces inside the spherical vessel of Fig. 1131 In spite of the wire connecting 
the inside and the outside spheres the field corrector introduced ensures that the electric field is spherically 
symmetric. 



